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We present numerical and analytical studies of a quantum computer proposed by the Yamamoto 
group in Phys. Rev. Lett. 89, 017901 (2002). The stable and quantum chaos regimes in the 
quantum computer hardware are identified as a function of magnetic field gradient and dipole-dipole 
couplings between qubits on a square lattice. It is shown that a strong magnetic field gradient leads 
to suppression of quantum chaos. 
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I. INTRODUCTION 

The fascinating idea of quantum computation stim- 
ulated numerous experimental efforts for realization of 
qubits based on various physical implementation of two- 
level quantum systems with controllable interactions (see 
e -g- H)- Liquid-state NMR implementation of quan- 
tum computation allowed to realize various quantum al- 
gorithms with few qubits including the Shor factoring al- 
gorithm and complex dynamics simulations Q . How- 
ever, the liquid-state NMR scheme does not allow in- 
dividual addressing of qubits and thus cannot lead to 
a large scale quantum computations The situation 
can be significantly improved in the case of solid-state 
NMR implementation of quantum computation Q . Cer- 
tain features of this scheme appear also in the all-silicon 
quantum computer recently proposed by the Yamamoto 
group 0, S H 111 ■ I n this proposal the qubits are spin- 
halves nuclei (e.g. isotopes 29 Si) placed on a 2D lattice 
on a surface of a crystalline solid matrix (e.g. of spin-0 
28 Si nuclei). A magnetic field gradient is assumed to be 
applied in the plane of the lattice to allow address qubits 
individually. Atpresent large gradients can be realized 
experimentally E3 and according to 0, || thousands 
of qubits can be addressed. In addition to qubit fre- 
quency gradient there are also dipole-dipole couplings be- 
tween qubits typical of liquid-state NMR 0, hj. These 
two elements, frequency gradient and dipole-dipole cou- 
plings between qubits, also appear in other proposals of 
quantum computers: e.g. for trapped polar molecules 
in electric field with gradient |l3j and trapped-ion spin 
molecules with magnetic gradient ^4|. Therefore the in- 
vestigation of generic properties of such systems is im- 
portant for future experimental implementations. 

Indeed, it is known that the residual couplings between 
qubits lead to emergence of quantum chaos and melting 
of quantum computer hardware |l5l Il6l Il7j| . It has been 
also shown [l8l Il9l | that these static imperfections give a 
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rapid decrease of fidelity of quantum computations and 
hence their analysis is important to preserve high com- 
putation accuracy. For a ID chain of qubits it has been 
shown that the introduction of frequency gradient gen- 
erally makes the quantum hardware more stable against 
static inter-qubit couplings. This result may have im- 
portant implications for proposals similar to those as in 
|l3llT3| . However, a more generic 2D case proposed in p| 
still requires a special study. Such a detailed study is the 
main aim of this paper. To analyze the generic properties 
the quantum computer hardware [j| we performed ex- 
tensive numerical simulations with up to 18 qubits. Our 
studies show that a sufficiently strong field gradient leads 
to a suppression of quantum chaos and emergence of in- 
tcgrablc regime where the real eigenstates become close 
to those of the ideal quantum computer without imper- 
fections. On the basis of obtained results we determine 
the border between this integrable regime and quantum 
chaos region with ergodic eigenstates for various system 
parameters. 

The structure of the paper is the following: in Section 
II we introduce typical physical parameters of the sys- 
tem, the results of numerical simulations and analytical 
estimates are presented in Section III and the discussion 
and conclusions are given in Section IV. 



II. DESCRIPTION OF QUANTUM COMPUTER 
MODEL 



Following the Yamamoto group proposal p| we in- 
troduce here a mathematical (YQC) model of a realis- 
tic quantum computer with field gradient and dipole- 
dipole couplings between qubits. For that we consider 
a two-dimensional array of n q nuclear spins 1/2 embed- 
ded in a crystalline solid, e.g. spin-1/2 31 P nuclei in a 
GaAs/AlGaAs quantum well or spin-1/2 29 Si nuclei in 
a matrix of spin-0 28 Si nuclei ||. In order to manipu- 
late the nuclear spins using radio-frequency (RF) fields 
each qubit should be distinguishable, i.e. a specific Lar- 
mor frequency has to be assigned to each of them. The 
Hamiltonian of this array of nuclear spins is the following 
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FIG. 1: Disordered spin lattice: Local deviation from the 
vertex position (n Xi b, n Vi b) is given by 5xi = 5 Xi b and 5yi — 
5 yi b (see text). The index i labels the spins and runs from 1 
to n q — n x n y . The magnetic field gradient has angle 9 with 
a;-axis (shown by arrow). 
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This Hamiltonian is widely used in NMR studies [Tlj . 



The operators /; 
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are the spin operators acting on the «th spin where a x ' v ' z 
are the usual Pauli matrices. The frequency luq + Suji is 
the Larmor frequency of the ith spin. Considering the 
spin- 1/2 nuclei as an ensemble of qubits, too + 8u>i is the 
energy between the two states |T) , |l) (or |0) , |1)) of a sin- 
gle qubit. The double sum runs over all the nuclei on the 
xy-plane and the term dij = d/qfj is the dipole coupling 
between the spins i and j. Here qtj = r^jh is the distance 
between two spins i and j in units of the spin lattice con- 
stant b; d is the coupling between two nearby spins on 
the square lattice (q^ = 1). In the Faraday geometry, i.e. 
the geometry where the magnetic field is orthogonal to 
the plane of the spin-1/2 nuclei (B = B z z), the coupling 
constant d is given by d = /x 7 2 ?i 2 /87r6 3 where 7 is the 
gyromagnetic ratio of the considered nuclei. A schematic 
view of the structure is shown in Fig. 1. In the present 
work we keep focus on the regime wq 3> StUi, d where the 
single qubit energy ujq is much bigger than the Larmor 
frequency shifts Suit and bigger than the dipole coupling 
dij between the qubits. This corresponds to typical ex- 
perimental conditions discussed in [3, Q ■ 

In a particular Larmor frequency shift Suji is as- 
signed to each nuclear spin. These Larmor frequency 
shifts are needed to distinguish one nuclear spin from 
the other, and can be seen as playing the same role as 
chemical shifts in liquid NMR Here, these Larmor 
frequency shifts originate from a magnetic field gradient 
g = V£>~ = gi lying on the xy-planc. For the ith nu- 



clear spin located at the position 
Larmor frequency shift is 
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where we define the characteristic frequency shift uj g = 
^bg. The Larmor frequency spacing between two nu- 
clear adjacent spins along the x (respectively the y) di- 
rection is tUg X = uj g cos 9 (respectively u> gy = uj g sin 9) . 
In Tab. 1 typical physical parameters are given for the 
case of 31 P and 29 Si nuclei. For example, if we consider 
a face-centered-cubic lattice of 31 P nuclei with nearest 
neighbor distance b = 3.9974A then d = 154Hz and 
the regime uj g /d = 10 in our simulations corresponds 
to u)g = 1540Hz which is equivalent to a magnetic field 
gradient g = 0.224T//im. If we consider now a lattice of 
29 Si nuclei with nearest neighbor distance b = IA then 
d = 2.4kHz and, for example, the regime (J g /d = 10 
in our simulations corresponds to u g = 24kHz which is 
equivalent to a magnetic field gradient g — 28T//im. 
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TABLE I: Typical physical parameters for different nuclei cor- 
responding to ui g /d = 10 in our simulations. 

For a case of square lattice wc consider the index i 
labeling row by row each nucleus in the array and thus 
running from 1 to n y (first row) then from n y + 1 to 2n y 
(second row) and this until i = n q = n x n y . The angle 9 
for the direction of magnetic field gradient is measured 
in radians and due to symmetry we assume 9 < 7r/4. Ex- 
perimentally, it is convenient to assign to the nuclei along 
the path indexed by i the Larmor frequencies in ascend- 
ing order. This ensures a total distinguishability of the 
nuclei. This condition is realized if n y oj g sm9 < uj g cos9. 
As the number of nuclei measured in a quantum dot can 
be as large 71rq — - Tidily r ~- / 10 4 , this condition becomes 
9 < riy 1 . Along the x direction the Larmor frequency 
spacing between nearest neighbors is then u gx ~ to g d 
and as a consequence the nuclei along that direction are 
highly distinguishable. The dipole coupling in is then 
negligible. Along the y direction the Larmor frequency 
spacing between nearest neighbors is u> gy ~ u g /n y <C tJ g , 
the nuclei along that direction are weakly distinguishable. 
With n y - 10 2 and using typical physical parameters of 
Tab. 1, we easily remark that along the y direction the 
dipole couplings dij cannot be neglected since d ~ w gy . 

In order to be more realistic we consider now the fact 
that the spins cannot form an ideal lattice, i.e. exper- 
imentally it is not possible to place each nucleus ex- 
actly on a regular rectangular lattice vertex. Consider 
the spin i at the vicinity of the lattice vertex (n Xi , n Vi ) . 
We define the deviations S Xi and S yi which character- 
ize the deviation of the spin i from the ideal vertex 
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position (n Xil n yi ). The position of the spin i is then 
r i = ( n xi + 5 Xi ) bx + (n yi + 5 yi ) by (see Fig. 1). We as- 
sume that the deviations 5 Xi and 5 Vi arc random and 
distributed in the interval [—8/2,6/2] where 5 is deter- 
mined by an unavoidable experimental error in the nu- 
clei positions. These spatial deviations modify weakly 
the dipole couplings dij between nuclear spins as it de- 
pends on the inverse of the third power of the spin-spin 
distance r»j . The main effect is the modification of the 
Larmor frequency shifts and as a consequence the Lar- 
mor frequency spacings between nuclear spins. For the 
ith spin the Larmor frequency shift is now 

Scot = (n Xi + 5 Xi )uj gx + (n Vi + 5 Vi )uj gy . (3) 

Thus the smallest Larmor frequency spacing between ad- 
jacent nuclear spins is 

uj' = min { 8u)i+ n — SuiA. . , . 

= (1-5) UJ gx - LO gy 5 

along the x direction and 

oj' gy = mm{6u i+1 -6w i } i=ltn9 _ 1 

= (1-5) UJgy - LOg X 5 

along the y direction. For 6 <G [0,7r/4] and for 5 < 0.5, 
u>' > which ensures that along the path i from one 
row to another the Larmor frequency increases. As 
we discussed before it is experimentally convenient that 
5u>i+i — 5u>i > for each couple of adjacent Larmor fre- 
quency inside a row, i.e. along the y direction. Us- 
ing lO this condition leads to the following inequality 
tan 8 > 5/(1 —5). Summarizing, for a large number n y 
of nuclei along the y direction and small spatial impre- 
cision i5 < 1 total distinguishability of the nuclei in the 
xy-plane will be achieved if the following condition is ful- 
filled 

5<6<n y l . (6) 

Note that the disorder introduced in J3J can also model 
the spatial inhomogeneity of the two-dimensional mag- 
netic field gradient g. Since the condition 5 < n y x has 
to be realized the condition of total distinguishability © 
imposes stringent experimental precision on the nuclei 
positions and on the two-dimensional magnetic field gra- 
dient. 

In this paper we consider a typical experimental situa- 
tion when the dipole couplings and the Larmor frequency 
shifts have comparable values d ~ lu 9 but ujq u> g ,d. 
In this regime, eigenstates of the Hamiltonian can 
be ordered by spin sectors (I z ) = Q^Zf). Eigenstates 
with the same (I z ) value form an energy band of width 
Sui ~ uj g y/n q . Within this band = n q \/(n q /2 + 

(I z ))\(n q /2 — (I z ))l eigenstates are contained. Nearest 
bands, i.e. bands with a difference of ±1 in (I z ), are well 
separated in energy since their spacing is ujq 3> 5lo. In 



this case, using spin flip operators = IJ ± ilj , the 
Hamiltonian Q can be rewritten as 

H = -ffdiag + -Hoff-diag (7) 

where 

n q 

tfdiag = £(^) + Sui) + 'I. j I: I',- (8) 

i=l i<j 

is the diagonal part of the Hamiltonian and where 

tfoff-diag =-\Y. ^ {# IT + ITI+) ( 9 ) 

i<j 

is the off-diagonal part. From © we straightforwardly 
remark that the off-diagonal part of the Hamiltonian (J7J 
is block diagonal, i.e. each block corresponds to a spin 
sector (I z ). This means that only noninteracting states 
with the same total spin value (I z ) are directly coupled 
by the dipole interaction. The diagonalization of any spin 
sector (I z ) can be performed independently of the others. 

The Hamiltonian J7J) determines the mathematical 
YQC model we study in this paper. Since the quantum 
chaos first appears in the middle of the spectrum |l5j | , we 
study the properties of the Hamiltonian (1) in the cen- 
tral band which is associated to the (I z ) = spin sector 
(respectively (I z ) = ±1/2 spin sector) for even (respec- 
tively odd) number of spins. This part of the spectrum is 
most sensitive to imperfections and emergence of quan- 
tum chaos first starts here. 



III. NUMERICAL RESULTS 

In order to characterize the emergence of quantum 
chaos in the YQC model J2J we use the level spac- 
ing distribution P(s) which clearly marks the transi- 
tion from integrable to ergodic eigenstates (see e.g. plj 
and Refs. therein). Here, s — (ei+i — ej)/A is the 
normalized level spacing for energy levels e, in a close 
vicinity of the middle of the central band (we use 10% 
of the spectral bandwidth 5oj); A = (e;+i — £j) is the 
mean level spacing in this vicinity. The form of P(s) 
is linked to the ergodic properties of eigenstates: an in- 
tegrable system has P(s) in the form of the Poissonian 
distribution Pp(s) = exp(— s), and P(s) is given by the 
Wigncr-Dyson distribution of the random matrix theory 
Pwd(s) = (tt/2) s exp(— 7rs 2 /4) if the eigenstates are er- 
godic. Fig. 2 shows typical examples of the level spacing 
distribution P(s) for a disorder strength 5 — 0.1 and for 
a magnetic field gradient angle 9 = 0.3rad. Each distri- 
bution P(s) in Fig. 2 have been calculated using 128700 
energy levels and Nd = 10 disorder configurations. As 
w g /d increases from u) g /d — to u) g /d = 10 the distribu- 
tion P(s) changes from the Wigner-Dyson distribution 
to the Poissonian distribution. Thus the YQC hardware 
is in the regime of quantum chaos for uj g < oj c g while for 
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FIG. 2: The level spacing distribution P(s) for a case with 
disorder fluctuations amplitude S = O.f, for 8 — 0.3, and 
for Lj g /d = 1 (•), 4 (■), 7 (□), and 10 (o). The solid 
curves show the Poissonian and the Wigner-Dyson distribu- 
tions. The number of disorder realizations is Nd = 10. There 
are n q — 16 qubits placed on the lattice of size 4x4. 



Lu g > uj g it is in the intcgrablc regime. According to the 
data of Fig. 2 the value of oj g is located in the interval 
4 < Lo c Jd< 7. 

To investigate the quantum chaos border uj g in more 
detail it is convenient to introduce the parameter |22j 



V = 



[P(s)- P WD (s)]ds 
Pwd(s)} ds 



(10) 



where sq ~ 0.4729 is the first intersection point between 
the two distributions Pp(s) and Pwd(s)- In this way 
r] tends to if the system is chaotic (P(s) — > Pwd(s)) 
and tends to 1 if the system is integrable (P(s) — > Pp(s)). 
This parameter ry is rather convenient for investigation of 
the transition from integrability to quantum chaos and it 
has been used in various studies of many-body quantum 
systems 

In Fig. 3 we show the dependence of 77 on the frequency 
shift ui g for different lattice sizes at a typical values of 
angle 9 = 0.3 and disorder strength S — 0.1. The tran- 
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FIG. 3: (a) Parameter rj as a function of uj g /d for different 
lattice sizes 5x2 (□), 6x2 (A), 7x2 (o), 4x3 (■), 5x3 
(•), 6 x 3 (*), and 4x4 (o). (b) Parameter 77 as a function 
sition from a regular regime (u g > u c g ) to a chaotic one G f the rescaled frequency shift uj g /uj c g where u) c g is defined by 
(u g < Ug) is clearly seen for each presented lattice size. n(oJg/d) = 0.3. For clarity only some curves presented in (a) 

are shown rescaled in (b). Upper inset: critical frequency shift 
uig as a function of n v , the straight line shows the linear fit (see 
text). Lower inset: critical frequency shift as a function of the 
total number of qubits n q = n x n y , where circles corresponds 



In order to monitor the transition more precisely by the 
dependence of uj c g on parameters, we define w g by the con- 
dition r](u)g/d) = 0.3. Fig. 3b shows 77 as a function of the 
rescaled frequency shift ajg/iOg. The lower inset of Fig. 3b 
clearly shows that uj c g is not a single- valued function of n q . 
This can be seen also in Fig. 3a where the lattice configu- 
rations 6x2 and 4x3 corresponding to the same number 



to n y = 2, triangles to n y = 3, and the square to n y — 4. 
In all cases 8 — 0.1 and 6 = 0.3. The number of disorder 
realization is Nd = 10 (Nd — 4 for the case 6x3). 



FIG. 4: (Color online) Entropy 5 n as a function of u) g /d, 
colors are proportional to the entropy. Red (gray) and blue 
(black) colors correspond to S n ~ 13 and S n = respec- 
tively; there are no eigenstates in the region with white color. 
Horizontal axis gives the value of u g /d from uo g /d = (left) 
to uj g /d — 20 (right). Vertical axis gives the energies of the 
central band eigenstates (arbitrary units). The size of the lat- 
tice is 4 x 4, n q — 16. Disorder corresponds to 5 = 0.1 and 
9 — 0.3rad. All the points on the figure has been computed 
for the same realization of the disorder. 

of qubits n q = 12 give two different curves of r\ and thus 
give two different values of uj g . The upper inset of Fig. 3b 
shows that ui g is satisfactory described by a linear depen- 
dence. A linear fit gives w g = \.2Qn y d. Also the main 
panel of Fig. 3b shows that the transition towards inte- 
grability becomes sharper as n y increases. Hence, for a 
100 x 100 array of 31 P spin-1/2 with a nearest-neighbour 
3ip_3ip distance b = 4A, a magnetic field gradient of at 
least g c = 2.8T//^m is needed in order to avoid emer- 
gence of quantum chaos in the static quantum computer 
hardware. 

The dependence found numerically 

u c g = Cdriy (11) 

with a constant C w 1.3 can be understand on the 
basis of the Aberg criterion |25j . This criterion has 
been proposed |25| to understand the conditions of emer- 
gence of quantum chaos and ergodicity in many-body 
quantum systems. The extensive numerical tests com- 
bined with a naly tical estimates have been performed in 
El IH El El IS E3 that allowed to establish 
the conditions of emergence of quantum chaos and dy- 
namical thermatization in variety of many-body quan- 
tum systems. According to these results the quantum 
chaos emerges when the coupling strength U becomes 



FIG. 5: (Color online) Entropy S n for uj g /d = 4 as a function 
of the magnetic field gradient angle 9. Colors are as in Fig. 4. 
Horizontal axis gives the value of 9 from 8 — (left) to 9 = 
7r/4 (right). Vertical axis gives the energies of the central band 
eigenstates (arbitrary units). Disorder strength is 8 — 0.1. 
All the points on the figure has been computed for the same 
realization of the disorder as in Fig4 at the lattice size 4x4. 



comparable to the spacing between directly coupled non- 
interacting states A c , i.e. when U ~ A c . It is important 
to stress that A c is exponentially larger than the level 
spacing between many-body quantum states which for 
YQC is of the order of A„ ~ uj g n q 2~ nq . 

In our model the noninteracting quantum register 
states are eigenstates \<f>i) of the diagonal Hamiltonian 
#dia g ©■ The off-diagonal part i?off-diag © contains 
flip-flop operators and couples these noninteracting states 
thus being responsible for the emergence of quantum 
chaos in the system. The ensemble of the noninteract- 
ing states {|</>i)} i=1 ff forms the quantum register ba- 
sis, i.e., the basis is composed of Njj = 2™' states with 
A^-dimensional vectors written in (|f) , ||)) or (|0) , |1)) 
representation (e.g. |^) = |0010100111... 01001)). The 
quantum register gives a convenient computational basis 
to perform quantum computations 0, each multiqubit 
state is then a linear superposition of quantum register 
vectors | 

To apply the Aberg criterion to the YQC model we 
note that the dipole interactions dij = d/qij vanish 
quickly with the interqubit distance and thus we can con- 
sider that only first (at most second) nearest neighbours 
qubits are coupled together. This means that U ~ d. 
For 8 <C 1 the states coupled by transitions i/ ff-diag © 
have a typical energy change lu 9 (assuming that only few 
rows in x contributes) and the number of such transitions 
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FIG. 6: Parameter 77 as a function of the magnetic field gra- 
dient angle 6 for uo g /d — 3 (•), 4 (■), and 7 (a). The disorder 
strength is 8 — 0.1. The lattice size is 5 x 3, Nd = 10. 

is of the order of n y . Thus A c ~ u g /n y that leads to the 
result . For 9 ~ 1 still the transitions in the direction 
perpendicular to the field gradient dominate so that in 
the relation (|11|) n v should be replaced by y/n^ that gives 

w c g = Cdy/n; ,(6~1) (12) 

with C ~ 1. 

The emergence of quantum chaos manifests itself not 
only in the level spacing statistics but also in the com- 
plexity of quantum eigenstatcs in presence of interac- 
tions. To demonstrate this fact in a quantitative way 
it is convenient to define the complexity of an eigen- 
vector \ip n ) by its entropy S n = ~Yln=i w ndog 2 w n i 
where w ni = 

K&HM ■ Here {\<j>i)} i=1 . N is the quan- 
tum register basis, or computational basis composed of 
Njj = 2 nq states. Here, iV^f-dimcnsional vectors are 
written in the (||) , |J.)) or (|0) , |1)) representation (e.g. 
\<j)i) = 10010100111... 01001)). The entropy S n quantifies 
the deviation of the eigenvector from a pure quan- 
tum register state. If S n = 0, \ip n ) = \4>n) is a quantum 
register state. If S n = n q , all the quantum register states 
\<pi) are equally present in the eigenvector \ip n }- As we 
particularly focus on the central band eigenstatcs, the 
maximum entropy attainable is S n = — log 2 (ij^/ 2 ]) < n q 
for an eigenstate \ip n ) equally composed by all the quan- 
tum registers | <j>f ) belonging to the (I z ) =0 (even case) 
or ±1 (odd case) sector. Indeed the symmetry of the 
Hamiltonian Q and the fact that loq 3> lo q > d guar- 
antee that the central band eigenstates \ip n ) can only be 
linear superposition of the | (f>f) quantum registers. 
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FIG. 7: Parameter r\ as a function of the disorder strength 
parameter 5 for ui g /d = 2 (•), 4 (■), and 6 (a). The magnetic 
field gradient angle is 8 — 0.3rad. The lattice size is 5 x 3, 
N d = 10. 



Fig. 4 shows S n as a function of (J g /d for a disorder 
strength S = 0.1 and a magnetic field gradient angle 
9 = 0.3rad. It is clearly seen that the entropy of eigen- 
states is reduced significantly with the increase of mag- 
netic field gradient. Thus a high gradient gives a suppres- 
sion of quantum chaos in the YQC model. The depen- 
dence of entropy on the angle 6 of magnetic field gradient 
is shown in Fig. 5. In agreement with the theoretical es- 
timates Ulljl. <|12f) there is no significant dependence on 
9. This is also true for the level spacing statistics as it is 
demonstrated in Fig. 6 where the parameter rj is clearly 
independent of 9. 

The dependence of r\ on the disorder strength d is pre- 
sented in Fig. 7. At weak and moderate values of rescalcd 
field gradient u> g /d an increase of S drives the system to 
a more integrable regime with higher values of r\. This 
tendency is similar to those found in [l5j j where disor- 
der also stabilized the integrable phase. However, in the 
YQC model this effect is weaker. Indeed, the increase of 
6 gives an increase of A c but not more than by a factor 
of two since there is always a contribution from a regu- 
lar field gradient. It is interesting to note that according 
to the data of Fig. 7 in absence of disorder the system 
is completely in the regime of quantum chaos (77 ks 0) 
at a small field gradient (ui g /d < 2). This means that 
at Lo g /d < 2 dipolc-dipolc interactions between qubits 
lead to a melting of quantum computer hardware and 
destruction of ideal computational basis. 

To determine the melting rate of the quantum com- 
puter hardware we compute numerically the local density 
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For strong couplings d (at T > A c ) the width starts to 
grow linearly with coupling r ~ driq 1 ^ similar to the 
case considered in [l^, [l|j, |23 • In this case pw (E) has a 
gaussian shape. The numerical data shown in Fig. 8 for 
18 qubits confirm these theoretical formulas. Indeed, for 
small ujg/d the width T is independent of uj g and pw{E) 
has a gaussian shape, while for larger u> g /d the distribu- 
tion pw(E) starts to have the Breit-Wigner shape and 
r oc l/ojg in agreement with l|14|) . For even larger val- 
ues of ujg/d one starts to enter in the integrable regime 
Lj g > u>g where the relation l|14|) is not valid anymore. 
The place of the transition is in a good agreement with 
the quantum chaos border l(TT|) which gives w c g = 3.9d 
(to be compared with the approximate value Lo c g f=a 5d 
obtained from data of Fig. 8). 

Thus an extensive amount of numerical data obtained 
with up to 18 qubits give good agreement with the the- 
oretical estimates derived for the YQC model. 



FIG. 8: Dependence of the width F of the local density of 
states pw on the frequency shift uo g for a disorder 5 = 0.1 
at a magnetic field gradient angle = 0.3rad and a number 
of qubit n q = n x n y = 5x3. The dashed line shows the 
gaussian regime with r oc d and the full line shows the Breit- 
Wigner regime with F oc d 2 /io g . Lower inset: an example of 
the local density of states pw for *-o g jd = 3; the full curve 
shows the best fit of the Breit-Wigner form with F — 2. Id. 
Upper inset: example of the local density of states pw for 
ujg/d — 1; the full curve shows the best Gaussian fit of width 
r = 2.52d. In both insets e/d is rescaled energy, crosses give 
numerical data. The transition towards integrability takes 
place near log (ui g /d) « 0.7 that corresponds to the quantum 
chaos border dill) . 

of states pw (see e.g. 0) defined as 

p w (E - Ei) =J2w ni S(E - E n ) (13) 

n 

where Ei is the energy of ideal quantum computational 
state (f>i given by ifdiag ©• The energy width T of 
this distribution gives the rate on emergence of quan- 
tum chaos 0, 0| . Thus in the regime of quantum chaos 
ujg < uig exponentially many states are mixed after a 
time scale 1/r (here H = 1). The number of mixed 
states is of the order of T/A n where A„ is a typical 
level spacing between many-body states. For weak cou- 
plings between qubits the distribution pw(E) has the 
Breit-Wigner shape with the width T ~ U 2 /A c (see e.g. 
PHH). For the YQC case this gives 

T ~ d 2 ^7T q /wg . (14) 



IV. CONCLUSION 

The numerical and analytical studies presented above 
establish the parameter regime where the ideal compu- 
tational basis of the quantum computer model (YQC) 
proposed in p| remains robust in respect to dipole-dipole 
couplings between qubits. They clearly show that a pres- 
ence of magnetic field gradient allows to suppress quan- 
tum chaos in the quantum computer hardware if the gra- 
dient g is larger than the quantum chaos border given by 
Eqs. I|ll[l . I|12|) . For typical parameters of Table I (e.g. 
29 Si, b = 1.9A) the YQC hardware is in the stable regime 
at g w 2 T//im for n q = 100 qubits and at g sa 20T/pm 
for n q = 10 4 qubits. These values can be realized with 
modern experimental methods. 

At the same time we should also note that the sta- 
bility of a quantum computer hardware is not necessary 
sufficient for high accuracy of quantum computations. 
Indeed, here we analyzed only the static properties of 
YQC. In a realistic YQC implementation it is also nec- 
essary to consider the accuracy of gate implementations 
and the effects of static imperfections on the accuracy of 
a concrete operating quantum algorithm (see e.g. [ljj). 
Future investigations are required to analyze such oper- 
ational accuracy of the YQC model p|. 

We thank T.D.Ladd and Y.Yamamoto for discussions 
of the specific properties of the quantum computer pro- 
posed in Q . This work was supported in part by the EC 
IST-FET project EDIQIP. 
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